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Abstract 

We complement a recent calculation (P.B. Gossiaux and the present 
authors, Ann. Phys. (N.Y.) in press) of the autocorrelation function 
of the conductance versus magnetic field strength for ballistic elec- 
tron transport through microstructures with the shape of a classically 
chaotic billiard coupled to ideal leads. The function depends on the 
total number M of channels and the parameter t which measures the 
difference in magnetic field strengths. We determine the leading terms 
in an asymptotic expansion for large t at fixed M, and for large M 
at fixed t/M. We compare our results and the ones obtained in the 
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previous paper with the squared Lorentzian suggested by semiclassical 
theory. 

PACS numbers: 72.20.My,05.45.+b,72.20.Dp 
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1 Introduction 



Lately, the transport of ballistic electrons through microstructures with the 
shape of a classically chaotic billiard has received much attention §. 
One observable investigated both experimentally and theoretically is the 
conductance autocorrelation function C(AB). With g(B) the dimension- 
less conductance, B the external magnetic field, and Sg(B) = g(B) — g(B) 
the difference between g(B) and its mean value, C is defined as C(AB) = 
5g(B)5g(B + AB). The bars indicate an average which experimentally is 
taken over the Fermi energy or the applied gate voltage and, in theories 
which simulate the chaotic billiard in terms of an ensemble of random matri- 
ces, over that ensemble. For ideal coupling between leads and microstructure, 
C depends only on the number M of channels in the leads, and on t, a mea- 
sure of AS defined below. 

The exact dependence of C on M and t is not known. In the semiclassical 
approximation, applicable for M>1,C was found to have the form of a 
squared Lorentzian ||. This result was confirmed by calculations using the 
supersymmetry technique j|, [|. However, the supersymmetric nonlinear 
sigma model leads to expressions for C(B) which so far have resisted all 
attempts at an exact evaluation. In the actual experiments pi, [7|, M is 
quite small. These facts and the desire to attain a deeper understanding of 
the supersymmetry technique prompted us || to calculate analytically the 
leading terms in the asymptotic expansion of C for small t at fixed M. In 
the present paper, we complement this study by calculating analytically the 
leading terms in the asymptotic expansion of C for large t and fixed M, and 
for large M and fixed t/M. 

Starting point is the random matrix model for the Hamiltonians of the 
chaotic microstructure for two values of the external field B. Our approach is 
identical to that taken in Ref . || . However, we introduce a new parametriza- 
tion of the coset manifold suitable for large values of t. We pay particular 
attention to the Efetov-Wegner terms (Refs. || to [[LlJ) generated by this 
procedure. 

In order to save space, we keep the introductory part as brief as possible 
and refer the reader to our previous paper 0. In Section |2|, we start with 
the form of the autocorrelation function obtained after averaging over the 
ensemble, after using the saddle-point approximation, and after integration 
over the massive modes. Unless otherwise stated, all quantities appearing 
in that section are defined as in Ref. M. The parametrization of the coset 
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manifold, and the integral theorem used in calculating the asymptotic for- 
mulae, are presented in Section |[ Simplifications of the integral theorem 
are discussed in Section |j. The asymptotic formula for large t at fixed M is 
derived in Section || Section ^| contains the conclusions. Additional mathe- 
matical details including the derivation of the asymptotic formula for large 
M at fixed t/M are given in several Appendices. 

2 Supersymmetric Form of the Autocorrela- 
tion Function 

The two magnetic field strengths are denoted by and B^ 2 \ We write 
g(B^) = with i = 1,2. The autocorrelation function 

C(t,M) = SgWSgW = ^V 2) - g^ g® (l) 

depends on M, the total number of channels in both leads, and on the pa- 
rameter t which is related to the area A of the billiard by \ft = k\B^ — 
B^\A/(2<p Q ). Here k is a numerical constant of order unity and <f)Q = hc/e 
is the elementary flux quantum |T7[. For sufficiently large values of B' 1 ' and 
we have = g^ = M/2, while g^g^ is given by || 

¥W ] = J V^T) exp < (Qr 3 ) 2 >) R{Q) detg~ M (l + QL). (2) 

Here, Q = T~ X LT. In the sequel, we use brackets < ... > as a short- 
hand notation for the graded trace. The matrices T belong to the coset 
space U(2, 2/4)/ U(2/2) x U(2/2) where U(2/2) x U(2/2) is the subgroup of 
matrices in U(2, 2/4) which commute with L. The symbol X>/z(T) denotes 
the invariant integration measure for the coset space. The 8N dimensional 
graded matrices L and r 3 are given by 

Larp,a'r'p' ( ) ^ S aa '8 rr iSpp' , (73) arp,a'r'p' ( ) ^ ^aa'^rr'^pp' • (3) 

The supersymmetry index a distinguishes ordinary complex (commuting) 
integration variables (a = 0) and Grassmannn (anticommuting) variables 
(a = 1). Later, the two values and 1 of the supersymmetry index a will 
also be denoted by b and /, respectively. The index r = 1, 2 refers to the two 
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fields. The index p = 1, 2 refers to the retarded and the advanced propagator, 
respectively. Finally, R(Q) denotes the source term given by the sum 

R(Q)=J2N j R j (Q) , (4) 

where Nt = M 2 /2, N 2 = M 4 /4, N 3 = iV 4 = M 3 /4, and where 

i?i(Q) = < GI™GI<& >< GlMGlW > , 
R 2 (Q) = < GI^GlM >< GI^GlW > , 

r 3 (q) = < gi^giWgiWgiM > , 

R 4 (Q) = < GIMGIMGIWGIW > , (5) 

with 

G(g) = (i + gL)- 1 . (6) 

The matrices 1^ are given by 

Ia'r'p',a"r"p" = ( — ^ 5 rr '5 pp '5 a ' a "5 r ' r "5 p ' p " . (7) 

The graded matrices Q can be parametrized in terms of 32 variables, half of 
them commuting, the others, anticommuting. In calculations of the average 
two-point function, one typically deals with a total of 16 integration vari- 
ables. Except for this increase in the number of variables and for the form 
of the source terms (which are, of course, specific to our problem), the form 
of our result in Eqs. (Jj) and @ is quite standard. In spite of this similarity, 
the increase in the number of variables renders a full analytical evaluation 
of the graded integral (0) very difficult. As pointed out in the Introduction, 
our analytical work is restricted to two asymptotic expansions of this inte- 
gral. We evaluate the leading term of the asymptotic expansion of C (t, M) 
in inverse powers of t and fixed M, and the leading term of the asymptotic 
expansion of C(t, M) in inverse powers of M and fixed t/M. Progress in this 
calculation depends crucially on the proper choice of the 32 variables used in 
parametrizing the matrices Q. 

3 Parametrization of the Saddle— Point 
Manifold 

In order to calculate the autocorrelation function for large t, we need a new 
parametrization of the saddle-point manifold different from the one used in 
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our previous paper H. Motivated by the paper by Altland, Iida and Efetov 



18], we write our coset matrices T as products of matrices Tj obtained by 



exponentiating the coset generators anticommuting with Ts, and matrices To 
obtained by exponentiating the coset generators commuting with T3, 



T = TjTo 



where T and Tj are the matrices 



(n 



11 



T 



,21 
'1 






f 12 
L l 
+ 22 
L l 






with 






'2 






f 12 
'-2 
22 



tf tf J 



Ti 







\tf 







j-22 
1 

f 12 
'1 



j21 
' i 

t 11 







22 



tf / 



(8) 



(9) 



t 



21 



for q = 1, 2, 3, 4 and 



tf = (l + tftff/\ 
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21 .12\l/2 
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In writing the matrices we label the rows of matrices of dimensions 2, 4, 
and 8 by the indices a, (pa) and (rpa), respectively. The indices follow in 
lexicographical order. The matrices of dimensions 4 and 8 are presented in 
block form. By construction, the matrix T (TO does (does not) commute 
with r 3 . Eq. (|9|) shows that the matrices T and Tj are constructed from the 
GUE coset matrices 



T, 



t zl t 
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22 



(12) 



the matrix To from the matrices T\ and T2, the matrix Tf from T3 and T4. The 
matrix elements {tf) aa i and their conjugates {tf)* aa i represent the Cartesian 
coordinates of T„. We denote this set of coordinates by x q . The invariant 



measure for integration over the coset matrices T = TjT has the form 

Vn(T) = l[Vfi G (T q ) V(T 3 ,T 4 ) , 



(13) 



where T>fiQ(T q ) is the invariant measure for integration over the GUE coset 



matrices T q , and P(T 3 ,T 4 ) is a function of eigenvalues of tftf and tftf 
given later. In Cartesian coordinates, 



vn G (T q (x q )) = nd(ti 2 ) aa ,d(t 



12 : 

g 1 



(14) 
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In spite of the simplicity of this expression, the Cartesian coordinates are not 
well suited for our calculation. Therefore, we follow Efetov || and change to 
polar coordinates. 



3.1 Polar Coordinates 

We introduce polar coordinates z r = (9%, (ft, , , jj, 7**) 7 2 , 7 2 *) f° r tne 
matrix To by the transformation 

ti 2 = uixi^uir 1 , tf = u 2 x\u i r y l (is) 

where u^, Xl 2 , X 21 are defined by 

u\ = exp ( ^ ^ , u 2 r = exp 7 J , (16) 

Xl 2 = i sin(0 r /2)e^ , X 21 = i sin(# r /2) e -^ , (17) 
6 r = diag(0 r ft , 91) , <p r = diag(#, g) . (18) 



The coordinates 0^,(f)^ are commuting variables, the cordinates 7^,7^* an- 
ticommuting variables. The relation t 21 = k(t^y implies X 21 = k(X) 2 )* so 
that 6 r = -k6* whereas (j) r = <j>*. Substituting Eqs. flT|) into Eqs. (|j), (|10|) 
yields 

T = UoAoU, 1 

where 



(19) 



U Q = diag(M 1 ,u 1 ,tt 2 ,M 2 ) , A 
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(20) 



with Xl 1 
and get 

where 



A 



22 



cos(9 r /2). We insert this expression for T in T = T/T 



T = U T I A U , 

T = Uq^TjUq 



(21) 
(22) 
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has the same form as the matrix Tj but with , t^ p replaced by 

tf = {ul)-Hf'4 , tf = . (23) 

The indices i are defined by i = 2, 1 for i = 1, 2 respectively. We replace 
the Cartesian coordinates x s of T s by the Cartesian coordinates x s of T s . 
The Berezinian of this coordinate transformation is equal to one. Since for 
s = 3, 4 the matrices t l2 i 2 ^ and t* 2 t 21 have the same eigenvalues, the measure 
has the same form in the old and in the new coordinates. We suppress the 
tildes. Finally we pass from the Cartesian coordinates x s of Tj to the polar 
coordinates z s = {6 b s , <fil,9f,<f)f,jl, 7**, 7 2 , 7 2 *) . We do this in the same way 
as in the transformation from x r to z r . Thus we get 

Tj = UjAjUf 1 , (24) 

where Ui,Ai denote the counterparts of Uq, A defined by 



/ M 1 



diag(«3,M4,M 



U 2 ) 



A, 



\3 








A 22 



A4 1 
Al 1 



V A 







X 12 



\22 
A 3 



(25) 



with A] 1 = A 22 = cos(# s /2). Collecting the results yields the matrix T as 
function T = T(z) of the polar coordinates z = (zi, z 2 , z 3 , 24), 



T = UoUjAjU^AoU^ 



-1 



(26) 



For all q, the matrices u v ql A* 2 and A 21 have the form shown in Eqs. ( JT6| ) and 
([T7|). By construction, the matrices U , A depend solely on the coordinates 
Z12 = (zi,z 2 ), the matrices Uj,Aj solely on the coordinates z u = (z 3 ,z 4 ). 
Moreover, Uq,Ui depend only on the anticommuting variables, A ,A/ only 
on the commuting variables. Substituting Eq. (^) into Q = T~ l LT and 
making use of the commutativity of Uq, Uj with L yields 



Q = UoA^UjAj'LAjUf'AoU, 



r-l 



(27) 



This is the expression of the matrix Q in polar coordinates. 
The integration measure dfi(z) = T>[i(T(z)) has the form 



dfi{z) = JJ„ d^c(^) p(8i 



3-1 ) 



(28) 
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where the d/i G (-2 g ) = 'Dp G {T q (z q )) denote the GUE measures given by 



dVG(Zq) = d[z q ]p G (9q) 



with 



d[ Zq ] = d[e q ]d[ Xq ] , PG (6 q ) = n« sin ^ < cos ^ >~ 2 

dftj = U a ^ , d|xJ = Il a <W R p dY q dY q * , 
and where p(6> 34 ) = V(T 3 (z 3 ), T 4 (z 4 )) denotes the density 

P(M = n a (^os^ + cos^) 2 n Q ^(cos^ Q + cos<)- 2 
/From these equations, 



dp(z) = d[z] U q PG(O q ) P(M , d[ Z ] = H q d[ Zq 



(29) 

(30) 
(31) 

(32) 

(33) 



The ordinary parts of the angles 9 q are integrated over the positive imaginary 
axis, the ordinary parts of the phases (J) q over an interval of length 2ir. The 
ordinary parts of Q{ are integrated over the interval (0,7r), those of d{ over 
the interval (0, 7r/2). 

An essential simplification of the calculation results from the observation 
that the matrices T possess a high degree of symmetry. We consider similarity 
transformations induced by the idempotent permutation matrices 
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(34) 



It is easy to see that Uq,A , Uj and Aj satisfy the symmetry relations 



Ji2U (z 12 )J 12 = U (w 2 i) , 
Ji2~Ui{z M )J 12 = U I (w 34 ) , 



Jl3^o(z 12 )Jl2 = A (u>2i) , 

Ji 2 A/(z 34 )Ji2 = A/(w 34 ) , 



(35) 



^34^0(^12)^34 = U (W 12 ) , ^34A (^12)^34 = A (W 12 ) , 

^34^/(^34)^34 = ^7(^43) , J 34 A 7 (z 34 ) J 34 = A/(w 43 ) , (36) 
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where w q = {6 b q) ip b , 0f, ipf, j3 q , j3 q , (3^, /3 2 ) represents a second set of polar co- 
ordinates defined by 

r q = -r q , P p q = t, P p q = ii- (37) 

For later developments it is important to consider the coordinates w q on 
the same footing as the coordinates z q . The Berezinian of the coordinate 
transformation from z to w is equal to unity, and the coordinates w are 
integrated with the same measure over the same domain as the cordinates 
z. Inserting the symmetry relations fl35p and (|36|) into Eq. (|26|) yields the 
symmetry relations 

T(z 12 ,Z U ) = Ji2T(w 2 i,w 34 )J 12 = ^34^(^12,^43)^34 • (38) 

Since the permutation matrices J12, J34 anticommute with the matrix L, the 
corresponding symmetry relations for the matrices Q = T~ l LT read 

Q(zn, Z34) = ~Jl2Q(W21, W U )J 12 = ~J?,aQ{w12, Wa 3 )J U • (39) 

3.2 Integral Theorem 

Eqs. (|29|) and (|30|) show that the measures dfic( z q) contain nonintegrable 
singularities located at < cos9 q >= 2(|A*| 2 + |A{| 2 ) = 2 < tftf >= 0. As 
we evaluate the integral 

X[F] = J Vfi(T)F(Q) (40) 

for some function F{Q) over the coset space in polar coordinates, the singu- 
larities cause the occurrence of additional terms (the Efetov-Wegner terms). 
These terms can be found by applying the method used in Ref. ||. We write 
the integral in Cartesian coordinates and exclude an infinitesimal neighbour- 
hood |(tg 2 )n| 2 + | C^- J 2 ) 22 1 2 < £ of singularities of dyUc^q) from the domain of 
integration. Making use of Berezin's theory of coordinate transforma- 
tions in superintegrals over functions which do not vanish on the boundary 
of the integration region, we change to polar coordinates, let e go to + , and 
obtain 

1[F\ = J d[z]p(e u ) H q (5(z q ) + PG (9 q )P(z q )) F(Q) . (41) 

Here S(z q ) denotes the S function of the polar coordinate z q , and P(z q ) is 
the projector on the subspace of functions of 4th order in the anticommuting 
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coordinates 7g,7g*,7g and 7^*. For a function f(z g ) of z q which is regu- 
lar at z q = and whose Taylor expansion in the anticommuting variables 
terminates with the 4th order term fi{z q ), we have 

d[z q ]5(z q ) f(z q ) = /(0) , P{z q ) f(z q ) = U{z q ) . (42) 



The product over the four factors (S(z q ) + pQ{9 q )P{z q ) s j in the integral (pLll) 
yields a sum of sixteen terms. We specify these terms by the indices a = 
(01020304). We put a q = (a q = 1) if the term contains S(z q ) (pG(O q )P(z q ), 
respectively). This yields 

1[F}= £ ?a m aAF} (43) 

ai,a2, £13,(24 =0,1 

where 

= / d[zK ia2a3a4 (tf 12 , m nf w p (^) f w) (44) 

with 

M=I1, Pg(WM • (45) 

The product denoted by (Ila ) extends over those values of q for which 
a q = (a q — 1, respectively) and equals unity if none of the a q s meets 
the definition. Eq. ( S3p is the desired decomposition of the integral ( PD| ) 



into surface and volume terms. The integral X im [F] extends over the entire 
coset space and represents the volume term. The remaining fifteen integrals 
represent the boundary or Efetov- Wegner terms. The term X ooo[-^] is equal 
to the value of F(Q(z)) at z = 0, Xoooo[-^] = F(L). For brevity, we use the 
symbol 2 a for the integrals X aia2a3(14 . 

4 The Efetov— Wegner Terms 

In our case, the integrand F in the integral 1[F\ has the form given in Eqs. (0) 
to ([jp. We write F as the product 

F = K RD , (46) 
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where K denotes the coupling term, 

K = e -m)<Qr,?> ; (47) 

R denotes the source term (f|), and D denotes the function 

D = detg- M (l + QL) = e -M<Hi+QL)> (4g) 

We show that for this form of F only four of the integrals F a [F] give a 
nonvanishing contribution. 

4.1 The Integrals X a [F] 

We rewrite Eq. (^7]) for Q in the form 

Q = UoA^QjAoU, 1 , Qj = UjAfLAjUf 1 (49) 

and use this form to rewrite the terms K, R, and D. For the coupling term 
K, we use that Uq,A and Uj commute with t%. As a consequence, we have 
< {Qr-i) 2 >=< (Qir 3 ) 2 >=< (A7 1 LA / r 3 ) 2 >. Using Eqs. © and Q, we 
find 

K = e -l<(QiTz?> = e 2tJ2 s <^e s > _ ( 5Q ) 

The exponent of the coupling term K thus depends only on the angles 9 S . 
The functions R and D both depend on Q via the matrix (1 + QL)^ 1 . We 
write Qi = L + SQi and use Eq. (^). We note that £/o commutes with L, 
introduce the matrix Go = (I + Aq 1 LAqL)~ 1 , note that Go commutes with 
Aq 1 and that GqAq 1 is diagonal and therefore commutes with L. We find 



where 
with 



1 + QLY 1 = U A ^0(1 + SW)- 1 AvUv 1 (51) 
SW = 5QiGqL (52) 



SQi = Qi-L, G = (1 + A^LAqL)- 1 . (53) 

Inserting the expression for (1 + QL)~ l in the defining Equation (|j) for i? we 
find that the anticommuting variables 7^, 7^* appear only via the matrices 

jtaO = c/-i/M[/ . (54) 
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The explicit form of these matrices as functions of 7*?, 7^* can be found with 
the help of Eqs. fllED and (|7|). We obtain 



nf *w nr p ^ PD = ic 5 ^ iC p (^) 5p 



where 



and 



(55) 
(56) 



Si 
S 2 
S3 
5 4 



with 



< VI^VI^VI^VI™ > , 



V 



(57) 
(58) 



The anticommuting variables 7^,7?* appear only in the matrix 5W . For D, 
we use Eqs. ©, (0), ©, © and © > and obtain 



where 



D = D e -M<Hi+8W)> ; 



-A/<ln(l+A ( 7 1 LA ( ,L)> _ -M £] r <ln(l+cos r )> 



(59) 



(60) 



depends only on 9 r . Again, the anticommuting variables appear only in 5W. 
For the integral I a [F] , we have 



(1) 



(61) 



with K, S and D given by Eqs. fl56|) and @. 

For later use we also give explicit expressions for SQi and Go m terms 
of# L 

m) yields 



,,. Yq and u'P. Using the formulae for U , A , Uj,Af in Eqs. (|20|) , (|25|) and 



(8Qf 



SQi 



5Q l 3 2 \ 

22 A/021 Q 



5Qi 2 5Q, 



(5gf5Q n 







21 4 ()' 5Q§ 2 y 



Go -2 



(G^-Gf \ 
Gf Gf 
G2 1 G2 2 
V GfGfJ 



(62) 
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where 5Q^ P \G^, P ' denote the matrices 

6Q] 1 = u\ (cos6 s - l)(ul)-\ 6Q] 2 = iu\ sin0 s e^(^)-\ 

SQf = -iu 2 s sin^e-^(^)- 1 , SQf =u 2 s (l- cos0 s )(^)-\ (63) 

Gf = itan(0 r /2)e^, Gf — ita,n(9 r /2)e~ i ' t "~, = 1 . (64) 

4.2 Symmetries of T a [F] 

The symmetry properties of our coset matrices imply symmetries of the inte- 
grals T a [F]. We use Eq. (|6l|) for X a [X], the symmetry properties in Eqs. (|35| ) 
of Uq, A , Ui, Aj, and we pass from the coordinates z to the coordinates w. 
We conclude that 

Sj(z 12 ,Zte) = Sjf(w2i,wu) (65) 
where f = 1, 2, 4, 3 for j = 1, 2, 3, 4 respectively. With Nj> = Nj , this yields 
S [z\ 2 , Z34) = 5 , (w 2 i,W34). Similarly, we conclude that 

D{z 12 ,z u ) = D(w 21 ,w u ) . (66) 

/From Eq. ((45|) we have p ai a 2 a 3 a 4 (0 12 , # 34 ) = p a2ai a 3 a 4 (0 21 ,9 u ), so that 



1a ia 2a 3 a 4 [F} = J d[w}p a2ai a 3 a 4 (9 2 i, 8 U ) 



x Ilf <*K) ilf p K) K(0u)S(w 21 , w u )D(w 21 , w u ) . (67) 

The coordinates W2 are integrated over the same domain as the coordi- 
nates W\, and the coordinates w q over the same domain as their counter- 
parts z q . Therefore, the integral on the r.h.s. of Eq. (^) equals the inte- 
gral l a2aia . ia4 [F]. Using the symmetry properties in Eqs. (^), we can re- 
peat the argument to show that S(zi 2 ,z 3 4) = S(wi 2 ,W4 3 ) and D(zi 2 ,z 34: ) = 
W43). We recall that p aia2 a 3 a 4 {di2, #34) — Pa^a^ifi 12 1 #43) and find 
that in this case the coordinate transformation yields the integral X aia2a4(l3 [X] . 
Combining these results shows that the integrals X Q [X] satisfy the symmetry 
relations 

X aia2a3a4 [X] X a2aia3a4 [X] X aia2a4a3 [X] X a2aia4a3 [X] . (68) 

These symmetry relations simplify the sum over the Efetov-Wegner terms in 
Eq. (H) to 

1[F] = X im [X] + 2X 1110 [X] + 2X 1011 [X] +X 1100 [X] 

+4X 10 io [X] + Xoon [X] + 2X 1000 [X] + 2Xooio [X] + X 000 o [X] . (69) 
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4.3 Vanishing X a [F] 



The sum in Eq. fl69|) simplifies even further: Due to the vanishing of some of 
the X a [F]'s, it reduces to four terms. We first consider the case a = (1110). 
The integration over 6(24) yields Z4 = 0. Using Eqs. ( |62"D to (|64|) for SQi and 
Go, we find that 



2 



/ 1 x x 

xi 1 

1 

V ulx-^uD^yi x 2 



u 3 X 3 {U 3 




X 2 
1 



and that 



detg(l + SW) = (1 - x{x{)(l - x b 3 x b 3) 
Here, x q , x q and y q denote the matrices 



x n 



itan(# g /2)e i09 , x q = itan(9 q /2)e 



-i<t>q 



Vq 



1 XgXq 



(70) 

(71) 
(72) 



The only part of ¥ which contributes to I^V Iq v ^ is the block V rPir y. 
According to Eq. ([TOD, ¥1,2,21 = ¥21,12 = 0, and the only nonzero part of S 
is the term N 2 S 2 . However, S 2 contains neither 73 nor 73*, and Eq. (|71"|) 
shows that the same is true for the function D. Thus P(z 3 )SD = 0, and the 
integral Z mo [-F] is equal to zero. For a = (1010), z 2 = z^ = and, therefore, 
¥12,21 = ¥21,12 = ¥21,22 = ¥22,21 = 0. Thus 5 = 0, and 2"ioio[-F] = 0. For 
the same reason, Ziooof-F 1 ] = 2ooio[-F] = ^oooof-F] = 0. The sum in Eq. 
simplifies to 



X[F] = J im [F] + 2T 1(m [F] + 1 1100 [F) + 1 0011 [F] 



(73) 



The term X 1100 [F] is independent of t and yields the limit t — > 00 of T[F]. 
For this term, K = 1, D = D , and S = N 2 S 2 = N 2 Y[ rp lrlr* U r < x r x r >■ 
Integrating over r ,7 r and 7* we are left with 

I noo [F] = (M 4 /4) J l[d[6 r } Y[ PG (6 r ) < tan 2 (# r /2) > D = M 2 /4 . (74) 



This is the disconnected part of g^gW. Thus, the autocorrelation function 
C(t, M) is given by 



C(t,M) —Z\\\\ [F] + ZZion [F] + loon [F] . 



(75) 
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5 The Limit of Large t at Fixed M 



In this section, we expand the autocorrelation function C (t, M) for fixed M 
asymptotically in inverse powers of t, and evaluate the leading term. In 
the limit f > 1, the integrals X a [F] appearing on the r.h.s. of Eq. (|75|) are 
dominated by the contribution of the neighbourhood of the surface 0% = 
where Qi = L, and where the graded trace in the exponent of the coupling 
term 

K = e -m<(Qr*r> ( 76 ) 

is equal to zero. 



5.1 Asymptotic Expansion 

We start from Eq. (|6lf) , introduce the rescaled angles 9" = i 1//2 #", and pass 
from the coordinates z s to the coordinates z s = {9 b s) cj) b s , 9[, 0{, 7* , 7^*, 7^, 7^*) 
with d[z s ] = dfz,]^ 1 and P(z a ) = P{z s ). We expand t~ 2 p a , K, V and D in 
powers of t~ 1//2 , 

t~ 2 Pa = r^ /2 p^ , K = t- nk/2 K ink) , 

VI (rp) = ^ t -n rp /2 V -(n rp )jCrp) ? £, = ^ t ~n d /2 D ^n d ) _ (77) 



n d 



The coefficients p^"\ K^ nk \ V^ nrp \ D^ nd) are functions of z r ,z s . Collecting 
terms of the same order in t^ 1 ^ 2 , we get 

UF] = Y, t ~ n/2l ^ ) i F ] ■ (78) 

n 

Here 1^ [F] denotes the sum of integrals 

jM[F] = ^r liw^-^lFj] (79) 

j n p n k nnni2n2in22n d 

with 

X (n P ^n d)[F]] = J H r d[ Zr )l[ s d[z s }pM 

X n! 0) 8(Zr) ^ ^) IL P{S a )K^NjS^D^ . (80) 
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Here, n v = (nun^i^)) and Sj 1 "' 1 are given by Eqs. fl5"T|) with Vio rp ^ 
replaced by V^I^ rp) , 



S[ nv) =< y( n n)j( n ">v( n ' 22 ' ) I^ 2) > < v( ni2) I^ 2) V( n2l) I^ 21) 



> 



S M =< y(n 11 ) I (n) v (n 12 ) I (l2) > < ^xJjCSDy^)^) > 
S^") =< ]/(mi)j( 1 l)y(n 2 2)j^22)y(n 2 i)j ( (21) 1/ („ 12 )j(12) > ^ 



The order of Sj is n a = J2 rp n rp- The sum over n p ,nk,n rp and n rf is 
restricted by the condition n p + rik + n a + rid = n. Since K is given by the 
exponential (|76|), all integrals contain the exponential factor 

K (o) = e 2E.<(^) 2 > . ( 8 2) 

On extending the domain of integration region over 6^ from zero to infinity, 
the series (f78|) yields an asymptotic expansion for X a [F] . Since p a , K, S and 
-D are even functions of 6f, only terms with even n appear, and the expansion 
proceeds in inverse powers of t. Thus, 

C(t, AO = E f J f m + 2Xffl [F] + zffi [F] ) . (83) 

n=l ^ ' 

To calculate the expansion coefficients Z d 2n \F{Q)\ , we use the symmetry 
properties of U , A , Uj, Aj. It follows (see Appendix |7.1|) that the integrals 
X^ l P' nk;nv;nd '[Fj] satisfy the symmetry relations 

n -(n p ;n k ;ni 1 n 1 2n2in22;n d )r j-r: <r( n P> n k;n22n2ini2nn;n < i) r j-, i 
• X aia 2 H FiJ - 1 a 2 a 1 ll [* j'\ 

^ r (n p ;n k -,ni2niin2 2 n2i;n d ) [ j-i i n -(n p ;n k ;n2in22niini2;n d ) r p n /o/f\ 
= 1 aia 2 ll l-fyj =i -a 2 aill [A?J> (°4J 

where j' = 1, 2, 4, 3 for j = 1, 2, 3, 4, respectively. 

To find the leading term of the expansion we recall two properties of 
Z(n P ;n k ;n v ;n d )[ F .y_ (j) Thig expression jg n0 nzero only when £ s P(z s )Sf v) D { - n ^ 

is nonzero and (ii) the expression does not vanish when integrated over the 
phases 0^. Inspecting the explicit expressions of the relevant low-order terms 
g( n v)jj(n d ) ag gj ven by the matrices (5Q pp 'Y n ^ and G pp ', we find that these 
conditions are first met for n a + n d = 4. Thus, the leading term C t> i(t, M) 
of the expansion is of order t~ 2 and given by 

C t>>1 (t,M) = t- 2 ( l[f n [F] + 21$ U [F) +4$n[F) ) • (85) 
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For each of the 2^ [F] in Eq. (§§), the sum on the r.h.s. of Eq. ( [75| ) reduces 
to the terms with n p = n& = 0. 



5.2 The Leading Term C t >i(t,M) 
We first calculate the volume integral 

Z*8i[*1 = E E ^rfr^^i] • (86) 

j ^11^12^21^22^ 

We suppress the tildes above z s and 9 S . Inspection of the explicit expressions 
of the relevant terms Sj D^ Ud ^ shows that the nonvanishing contributions 

to J2 S P{z s )S ( j lv " > stem solely from the terms which are linear in all four 
matrices (8Q l2 )^\ (SQ 21 )^ and which for each r contain the same number of 
G) 2 and of G 21 . We restrict the sum over X^if"" [Fj\ accordingly, employ 
the symmetry relations Eq. fl84|), and find 



T ( 4 ) \m - T (0;0;llll;0)r pl 
■illllFj - -^1111 

+J (0,0;0000;4) ^ + ^OgOOO^ ^ + 2J (0,0;0022;0) ^ 

+2J (0;0;0220;0) ^ + 2 ^OgOlOl^) ^ + 4J (0;0;0121;0) ^ ( g ?) 

The integrals X^n^ [-P}] have the form shown in Eq. (j30|). Explicit expres- 
sions in terms of polar coordinates for the V^ n \ D^ n > contributing to Eq. fl8"7p 
are given in Apendix [7.2| . Substituting these formulae and applying the pro- 
jectors P(z r ) yields 

4t\i[F] = (M 2 /64) J dizjp^K^ Dollar P 

x( 2 < xPy 2 X^ yi >< X i i ) y 1 X^y 2 > 
+M 4 < Vl >< y 2 >< X^^X^ca >< x£ ) x 1 Xl 1) X2 > 
+4M 3 < yi >< xP^X^xi >< x£ ) x 1 xi 1) X2 > 
-2M 2 < yi >< X ( i ) y 2 X ( i ) x 1 X ( i ) y 2 X { l ) x 1 > 
+2M 2 < X^x 2 xi 1) x 1 y 1 >< X^x 1 xi 1) x 2 y2 > 
+2M 2 < X^x 2 y 2 xl 1) x 1 y 1 >< x£ ) x 1 X<l ) x 2 > 

-AM < X^y^x^y^xm > ) , (81 
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where Xf^,Xf' denote the matrices 

x^ = {5Q^r = lu \e s ^(ul)-\ 

X™ = ( 5 Q? )W = -in 2 e^iul)- 1 , (89) 

and where 

are the matrices introduced in Eq. fl72|). The density Pirn is given by 

pfAi = nPG^)iiPG> s ) , p^m = < m 2 >- 2 , (9i) 

r s a 

and P denotes the projector P = Y[ s P( z s) which we suppress in the sequel. 
First we integrate over </> s ,7 s ,7*. This can be done making use of 

/ 1] d[ Xs ] < X« Aj X^ Bj>< Xi 1] A H xf> B u > 

J s 

= II < (°s) 2 ><A r >< B I >< A n >< B„ > , 

s 

f n &[Xs] < XP Aj Xi 1] Bj>< XP A H Jtf B u > 

J s 

= n < (^) 2 > < A i B n >< B i A n > > 

s 

f n d[ Xs ] < xg> Aj X« Bj A H x| 1} B n > 

J s 

= II < i 9 *) 2 >< A i >< A n >< BjBu > . (92) 

s 

These formulae are valid for all two-dimensional graded matrices Aj,Bi, 
An, B n which are independent of 4> s ,ls,ls^ an d can be verified with the 
help of 

/ d^rwu«M)-u = R p+a2 ww^r 1 • (93) 

Integrating over <f) r ,lr, 7* then simplifies the volume term to an "eigenvalue" 
integral over the angles 6® , 

4UF] = (M 2 /64) / II A dA PnnK^D U < (6 S ) 2 > 



x[< yi > 2 ((M 4 + 2M 2 + 2) < y 2 > 2 -AM 3 < y 2 (l - y 2 ) > 

-4M< Vl (l -y l ) > ( < y 2 > 2 -M< y 2 {l -y 2 )>)y (94) 
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Since the exponentials , D factorize into products of exponentials each 
depending on one of the 9 q , the integration over angles simplifies to the 
evaluation of two-dimensional integrals over 6 q . The integrals over 6 r can be 
done as in Ref . . With 

J dK]p G (^)e- M<ln(1+coser)> < y r > 2 = (1 - M 2 )" 1 , 

J dK] PG (^)e- M<ln ( 1+cos ^) > < y r (l -y r )> = M(l - M 2 )" 1 (95) 



and 



P % OO P OO 

= 4 / d9 b 3 / d^e 2<(es)2> ^f < {6 S ) 2 >~ l = 1/2 , (96) 
Jo Jo 



the result is 



K4) r ^ M 2 M 4 -2M 2 + 2 



Jllll[F] - 256 (M*- 1)2 • (97) 

The boundary terms can be evaluated in the same way, see Appendix |7.2| . 
For a = (1011) we have Z2 = 0, which reduces the number of contributing 
integrals. We find 

T ( 4 ) r PI - t(° ;0;1111; °) [PI i T (0;0;0022;0) r p , 9T (0;0;0121;0) r p , 

_M 2 M 2 -2 

"256M^T • (98) 
For a = (0011) we have = 22 = 0, and only one integral contributes, 

T ( 4 ) rpl _ r (0;0;llll;0)r pl _ M 2 

-Ami i-n - Axni ^ ^ ~ 128 ' ^ ' 

We thus find that the leading term is 

C t>>1 (t , M) = t- 2 ( J^n [F] + 2J{ 4 ) 11 [F] + jg^ [F] ) 
1 M 4 /M\ 2 



- 



4(M 2 -l) 2 V8t/ ' 
This expression is discussed in the next Section |[ 



(100) 



20 



6 Summary and Conclusions 



We have investigated the magneto-conductance autocorrelation function for 
ballistic electron transport through microstructures having the form of a 
classically chaotic billiard. The structures were assumed to be connected to 
ideal leads carrying few channels. Assuming ideal coupling between leads and 
billiard, we have described the system in terms of a random matrix model. 

The autocorrelation function depends only on the field parameter t, spec- 
ified by the field (flux) difference, and on the channel number M. Using the 
supersymmetry technique, we have calculated analytically the leading terms 
in the asymptotic expansion for large t at fixed M, and, in Appendix |7.3.2|, for 



large M at fixed t/M. To this end, we have employed a new parametrization 
of the coset space which is particularly convenient for the present problem of 
partly broken symmetry. Applying Berezin's theory of coordinate transfor- 
mations in superintegrals, we succeeded in formulating the relevant integral 
theorem, and in identifying volume and boundary (Efetov-Wegner) terms. 
Symmetry properties of the coset matrices helped us in considerably simpli- 
fying the evaluation of these terms. We believe that the method developed 
in this paper is of general interest for the supersymmetry technique, and we 
hope it will be helpful in other cases. We have shown that the Efetov-Wegner 
terms are essential for details of the large t behaviour at small M. 

To discuss our results, we construct a function C$ L (t,M) of the form 
of a squared Lorentzian which has the same leading term at large t as our 
result (|100P . We use the value of the autocorrelation function at the origin 
t = given by the variance M 2 /[4(M 2 — 1)] of the dimensionless conductance 
g(B), see Ref. This yields 




c L(t,M) = ———[ i + —VW^T ) . (loi) 



In our earlier paper ||, we found that the leading terms in the asymptotic 
expansion of the autocorrelation function for small t are 

M 2 4M 3 
C[t ' M) ~ 4(M 2 - 1) " (M 2 - l) 2 * 

16M 2 (3M 4 - 11M 2 + 36) 2 3 
+ (M 2 - l) 2 (M 2 - 4) (M 2 - 9) 1 ( 102 > 
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The squared Lorentzian consistent with the first two terms is 

C-- L (t,M)= ,f Jl + EJf-)- 1 . (103) 
SLK ' ; 4(M 2 — 1) \ M M 2 — 1 / V ; 

The two squared Lorentzians are different. This shows that we cannot find a 
function with the shape of a squared Lorentzian which has the same leading 
terms as the exact correlation function for both small and large values of 
t. In other words, the exact shape of the autocorrelation function is not a 
squared Lorentzian. However, with increasing M, the differences between 
Cg L (t, M) and C|£(t, M) decrease, and for M ^> 1, the semiclassical result 

C ^(*> M ) = i( 1 + ^) (104) 

does have the leading terms required by Eqs. ( |100| ) and ( |102j ) . Moreover, 
the numerical differences between the three functions Csl,Csl an d C|£ are 
small even at small M. ^From the numerical simulations discussed in Ref. M 
the difference between the exact RMT autocorrelation function and the best 
squared Lorentzian fit does not exceed 5% in magnitude even for the smallest 
value of M, M = 2. 
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7 Appendices 

7.1 Symmetry Properties of the Integrals 

^-(n p ;n k ;n v ;n d ) r ^ -, 
Ana 2 ll [^iJ 

Similar to the symmetry properties of I a [F] discussed in Subsection [4.2|, these 



properties follow from the symmetries of the matrices Uq, A , Uj, A/ under the 
similarity transformations by the permutation matrices J 12 and J34, and can 
be found in the same way. The integrals are given by Eq. (|8D|). We use 
Eqs. (|35|) and change from the coordinates z r ,z s to the coordinates w r ,w s , 
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where w s = {6 b r1 ip b Tl 6l ,ip£ , Pi, Pi, Pi, Pi) • Since the expansion terms satisfy 
the same relations as the expanded functions, we have 



ri(niini2n2in22)/'^ ~ \ 



q[7 



(n22".2ini2nn) 



w 2 i,w M 



with f = 1, 2, 4, 3 for j = 1, 2, 3, 4, respectively. Similarly, 

D^(z 12 ,z u )=D^\w 21 ,w u ) . 



(105) 



(106) 



With p££n(0l 2 ,034) 



^20111(^21,^34) we finally find that 



(n p ;n k ;nnni2n2iri22;n d ) 



aia2ll 



ft-: 



n dH n « p£ani(^i, m n! 0) w 



(ri22n2ini2nii) 



^i^^W,^). (107) 



Since the coordinates w 2 are integrated over the same domain as W\, and the 
coordinates w r ,w s over the same domain the integral on the r.h.s. 



of Eq. QIU7D is the integral Z££ : 



(n p ;n k ;n22n2in 1 2n 11 ;n d ) 
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[Fji] . The same procedure 



can be repeated using the symmetry properties with respect to J34 instead 
of those with respect to J\ 2 . In this case, the coordinate transformation 
yields the integral x^^ Ul2niin22n21,ni \Fji]. Summarizing and combining 
the results we find that the integrals F^'^ k l ' nr ' nd \Fj} satisfy the symmetry 
relations in Eqs. (jS4j). 



7.2 The Integrals 1^[F] 

For the volume termX^jF], the functions V^ 2 \ D^ 2 \ take the form 

/ 



v<n = - 

4 

























-x? 





x? yi 





x? Vl 








x? 


x 2 X ( l ) y 1 











X?y 2 \ 



/ 



X^x 1 xf ) y 2 






A 4 X1A3 y 2 



23 



= (MD^/A) < Xi 1) x l xi 1) x 2 > , 
= (M^^/ie) < X^xzXPx! >< xi 1) x 1 X ( l ) x 2 > , 

where Xjj 1 ', Xf' are the matrices given in Eq. (|89D and x r ,x r , y r the matrices 
introduced in Eq. ([721). For brevity, we present here only those parts of 
£)( 2 ) j £)( 4 ) which give a nonvanishing contribution to the integral. 

For the boundary terms, inserting the matrices V^ n ' and applying the 
projectors n£ P( z r) yields 

1%[F] = (M"/ 256) J d^mA^nK^Doll^ 

x ( -2 < xPkXPyx >< X^yxX^k > 

+ M 2 < yi >< X^kX^xxX^kX^xx > 

+ 2M < X^kX^xxX^kX^yxXx > ) , 

l&i[F] = (M 2 /512) J Y[ s d[z s }AxK^D 

x < XpkX^k >< Xl 1] kxi 1] k > , (108) 

with p(°) = Pc(^r) HspG^i^s) ■ The phases and the anticommuting vari- 
ables can be integrated out using the integration formulae (^). This simpli- 
fies the boundary terms to the eigenvalue integrals 

lx%[F] = (M 2 /64) J dftllLdpJ AiK^Dol[ s < [9 s f > 
x ( _ (m 2 + 2) < yi > 2 +2M < - j/i) > ) , 
= (M 2 /32) J UA 9 ^ Ki0) H < > ■ ( 109 ) 
These integrals can be done using Eqs. and 



7.3 The Limit of Large M at Fixed t/M 

This limit has been considered, and the leading term worked out, previously 
in Refs. ||, [5|. Thus, the present Appendix provides a test of our proce- 
dure. At the same time, comparison between the two approaches shows the 
relative simplicity of the present approach. We proceed in the same way as 
in Section |5|. Many steps are similar and will not be repeated. We expand 
the autocorrelation function C(t,M) in inverse powers of M and evaluate 
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the leading term. In the limit where M, t ^> 1 at t/M fixed and arbi- 
trary, the integrals X a [F] appearing in the sum ([F^), i.e. the integrals with 
a = (1111), (1011), (0011), are dominated by the contribution of the neigh- 
bourhood of the surface 9" = 9® = where Q = L and where the graded 
traces in the exponents of 

K = e -(t/2)<(Qr 3 ) 2 > ^ jj = e -M<ln(l+QL)> /-qq) 

are equal to zero. 

7.3.1 Asymptotic Expansion 

We start from Eq. (|6lD , rescale the angles 9 q = M 1 ^ 2 9 q i , and pass from the 

coordinates z q to the coordinates z q = (9 q ,<p q ,9^,<p^, r y q ,j q *,j q , r y q *) , with 

d[z}p a U^ S(z r ) = d[z]M~ l p a nl 0) 5(z r ) for I = E g a q . At fixed t/M, we 
expand M~ l p a , K, V and D in powers of M -1 / 2 . The expansion coefficients 
p{n P ) ^(n k ) ^y(n rp ) ^ ]j{n d ) are f unc ti ons Q f % q . We take into account that the 

Nj also depend on M, and write Nj = M U ^ 2 N^ 3 \ with Vj = 4,8,6,6 for 
j = 1,2,3,4, respectively. We collect all terms of the same order in M -1 / 2 
and get 

la{F] = Y,M~ n/2 li n) {F} . (Ill) 

n 

Here X^> [F] denotes the sum of integrals 

l^ n) [F]=Yl Yl j{n P ; nfc in„ W )^.] ( m ) 

j n p n k n 11 n 1 2n2iri22n d 

with 

X {n p ;n k ;n v ;n d :u j )^ = J d [ 5 ] p K) 

x Ilf *&) 11? P &) U s P{*.) K^N^S^D^ . (113) 

For fixed j, the sum over n p ,nk,n rp ,nd is restricted by the condition n p + 
nk + n a + rid = n + Vj- All integrals contain the exponential factor K^D^ 
with 

K (0) = e (2t/M)J2 s <(8s) 2 > ^ £,(0) = e (V4)E 9 <(^) 2 > . (H4) 

On extending the domain of integration region over 9? from zero to infinity, 
the series ( |111| ) yields an asymptotic expansion of 1 a [F] . Since p a , K, S 
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and D are even functions of 9®, only the terms with even n appear, and the 
expansion proceeds in inverse powers of M, 

OO , s 

C(t, M) = Y, M ~ n J mi ^ + 2J i ( on [F] + [F] • (1 15) 

n=l ^ ' 

The integrals if"^^-"^^^] sat i s f y the same symmetry relations as the 
integrals X^ lp ' ,nk ' nv ' nd ' ) [Fj], cf. Eq. (f34j). The leading term of the expansion is 
of the order M° , and is given by the sum 

C M>>1 (t,M) = l$ n [F] + + Ji° ) 11 [F] . (116) 

7.3.2 The Leading Term C M >i(*,M) 

The nonvanishing projections J2 S P{z s )S^ av ^ stem solely from the terms 
which are linear in all four matrices (SQl 2 )^\ (SQ 21 )^ and which for each 
r contain the same number of (Cr* 2 )^ and of (G 21 )^. We employ the sym- 
metry properties of x { a p ' nk ' nv ' nd ^ ) [Fj\ and find that 

T (0) rpl _ T (0;0;llll;4:8)r pl 

•^miL^J - -Mill [*2\ 

+ 4J^ m3;2:8) [F 2 ] + 2I^ 31 ^ [F 2 ] + 2Xff° ;1111;2:6) [Fs] , (117) 



T (o) 



ion[*1 

(0) 



X 



^0011 — X, 



(0;0;1133;0:8) 
1011 

(0;0;1111;0:4) 
0011 



[Fi] + 21^1 
[Fi). 



(0;0;1131;0:6) 



[F 3 



(118) 



For the volume term [F] , 



o 

2xS 1} 




2x« 



(i) 



(i) 



-X, 



2s« 



4M 



2x 



(i) 



\Z(3) = — 


















y(l) T (l)y(l) 
A3 A 4 







A 4 A 3 
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D® = p(°> /4) < xPx^X^x? > , 
D W = p(o)/i 6 ) < X^x^X^xP >< X^^X?^? > , 

where X^, X,W are the matrices introduced in Eq. (159), and I > denote 
the matrices 

4 1 ) = (G T 12 ) (1) = i{6 r /2)e^, x« = (G 21 ) (1) = z(^/2)e"^ . (119) 

We present only those parts of V^,Z)^,£)^ which give nonzero contri- 
bution to the integrals. The densities p^ have for all a the form p£) = 
rig 1 ^ Pg (^g)j whh Pq* given in Eq. fl9T|). The integrals are evaluated in the 
same way as in Subsection p^2| . The integration over the eigenvalues is done 
making use of 

/ d[e r ) P <% (e r )e^ 2 > < y ?) > = i , 

/ d[^]p ( 2 (^)ei<^) 2 > < y( 2 > > 2 = -1 (120) 

and 

d[e s ]p {0 i (e s )e^ 1+ W M V< e ^ 2> < (6 S ) 2 > = 4[l+ (%t/M)Y l . (121) 



The contributions stemming from the integrals appearing in the second row 
of Eq. ( |117[ ) cancel mutually, as well as the contributions of the two boundary 
terms. Thus the evaluation of the leading term reduces to the evaluation of 
the integral 



xST nm [F 2 ]=jMp^ n K^ 

x X 2 (8) < /C") V V 4 12 ) >< V ^ 4 21) V™ 4 22) > 
= (1/64) / d^p^K^D^ n r < *?> > 2 



x 



IL^ < jtf^W >< xPxPxPxP > 



= (V64) / n, d^pf^D^ n r < v? > 2 il < ^? > ■ 

The integrals over 6 r yield the unit diffuson factor independent of t/M, the 
integrals over 6 S yield the cooperon factor 16^1 + (8t/M)^j . The result is 

CM>>i(t,M) = ~(l + ^) . (123) 

27 



This result agrees with the one obtained by Efetov and by Frahm || in 
the framework of supersymmetry, and with the squared Lorentzian found in 
the semiclassical approach ||. 
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